Physical limits on computation by assemblies of allosteric proteins 
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Assemblies of allosteric proteins, nano-scale Brownian computers, are the principle information 
processing devices in biology. The troponin C-troponin I (TnC-Tnl) complex, the Ca 2+ -sensitive 
regulatory switch of the heart, is a paradigm for Brownian computation. TnC and Tnl specialize 
in sensing (reading) and reporting (writing) tasks of computation. We have examined this complex 
using a newly developed phenomenological model of allostery. Nearest-neighbor-limited interactions 
among members of the assembly place previously unrecognized constrains the topology of the sys- 
tem's free energy landscape and generate degenerate transition probabilities. As a result, signaling 
fidelity and deactivation kinetics can not be simultaneously optimized. This trade-off places an 
upper limit on the rate of information processing by assemblies of allosteric proteins that couple to 
a single ligand chemical bath. 

PACS numbers: 87.16.Xa, 87.15. km, 02.50.Ga, 89.75.-k 

Keywords: allosteric regulation, signal transduction, energy landscape, Markov network, Markov chain, 
statistical mechanics 



Introduction Cells continuously regulate their function by responding to changing concentrations of diffusing 
molecules, called ligands. Ligands are sensed by allosteric proteins — measuring devices that communicate ligand 
binding to an input domain as a change in structure [TJ |5] or dynamics [5] in a distinct output domain. Often, 
allosteric proteins arc organized into allosteric supramolecular assemblies 4, 5 , where function is modularized with 
proteins specializing in ligand sensing or reporting. These assemblies operate as digital logic buffers, with a binary 
input (bound- or unbound-ligand to the input of the sensor) and binary output (bi-metastable state [B] of the reporter 
protein). Important examples include the family of G-protein coupled receptors and the NF-kB family of transcrip- 
tion factors [7] ■ Brownian computers are a class of computational models that use random collisions to stochastically 
explore the low energy portions of the computer's configuration space [5]. The allosteric supramolecular assembly is 
a paradigm for Brownian computation. 

Inquiry into the physical limits of computation began with the Maxwell demon (MD), a thought experiment by 
James Maxwell [9]. The MD was the first model to include both tasks of computation: sensing and reporting, also 
called reading and writing. Szilard captured the function of the MD in a computationally tractable model of an 
engine operating on a one dimensional gas [10 . The Szillard engine drew attention to two physical concepts in 
computation: the relationship between information and entropy and the problem of system resetting. Both sensing 
and reporting require a physical memory. Landauer, representing condensed phase memory as a symmetric bi-stable 
potential, argued that system resetting has a minimum thermodynamic (entropic) cost of fee In 2 dissipated per bit 
processed — Landauer's principle [8j[TT]. The bi-stable potential model applies to a restricted class of simple systems, 
which limits the ability of Landauer's principle to reveal the physical limits of computation [T5]. In particular, the 
bi-stable potential does not adequately model the complex architecture needed to perform the sensing and reporting 
tasks of computation. Other attempts to place physical limits on computation |13j recognize the role of kinetics but 
ignore the issue of system resetting. 

Here, we explore how system complexity and system resetting impose physical limits on the rate of computation by 
assemblies of allosteric proteins. Our findings are based on the free energy landscape of the two-component cardiac 
regulatory assembly, a Brownian logic buffer. The cardiac regulatory assembly is a Ca 2+ -sensitive switch that allows 
the heart to undergo periodic contraction and relaxation necessary for its pumping action [14] . The switch consists of 
the sensor protein troponin C (TnC), a receptor for Ca 2+ , and the reporter protein troponin I (Tnl), a regulator of 
muscle contraction. Function is described in terms of the system's free energy landscape, resolved at the mesoscopic 
level. The landscape must support both activation and deactivation (resetting) components of the signaling cycle. 
The topology of the landscape is constrained when communication between the system components — Ca 2+ , sensor, 
and reporter — is limited to nearest-neighbor interactions. We show that in the constrainted landscape, increases in 
signaling fidelity are offset by decreased resetting speed. This trade-off between signaling fidelity and resetting speed 
can limit the rate of computation by Brownian computers. 

Signaling Dynamics — The essential statistics and dynamics of the cardiac regulatory assembly are described by 
a mesoscopic phenomenological model of allostery. System dynamics involve stochastic jumps between metastable 
configurations of the assembly, called system-states. Each system-state Si is represented by a unique three bit binary 
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string (soW s 2)- Bit s gives the liganded status of the sensor (s : = unbound, 1 = bound), si is the output state 
of the sensor, and S2 is the output state of the reporter (sj>o: = inactive, 1 = active). The output state of the 
assembly is determined by the state of the reporter s 2 . The configurational phase space, ft = {Sj : i = 1, 2, . . . , 8} = 
{(0\00), (0\10), . . . , (1\11)}, consists of eight system-states. Only single bit Sj transitions are allowed. Transitions 
are represented by <7+ = (sj : — * 1), and aj = (sj : 1 — ► 0). In discrete time increments, r = tvo, where vq is 
the fastest barrier crossing attempt frequency of all system transitions (the clock frequency) , the system evolves as a 
Markov chain 

n 

P J -(r + l)=£;n ii (G(/ i (r) J /OP i (r), (1) 
»=i 

subject to the initial distribution {Pj(0)}. P»(r) is the probability that the system is in state Sj at time t. System-state 
transitions are governed by transition probabilities IX^ = Ii(Sj <— Sj) that depend on the free energy landscape of the 
system G(/i(r)), which is a function of the time-dependent chemical potential of the Ca 2+ ligand /i = fcsTln[Ca 2+ ] 
and time-independent friction coefficients /j for the system-states Sj [2 [16]. is relative to the standard state, 
taken as 1 M [Ca 2+ ]. The free energy landscape of the system is the collection of free energy surfaces along each 
elementary transition. The Markov network is a graph associated with where the Sj are represented as vertices 
and elementary transitions are represented as edges. 

Fig. [T^ shows the basin-limited free energy landscape — a simplified representation of the free energy landscape — of 
a generic two-component allosteric assembly. The free energies of the system-states {G(Sj)} are shown on the z-axis of 
the Markov network with transitions in S2, Si and sq along the x, y and z axes. By partitioning the {G(Sj)} into the 
ligand-unbound {G(0\siS2)} and ligand-bound {G(l\siS2)} surfaces, the basin-limited landscape provides a visual 
description of how the system's energetics supports signaling. Activation and deactivation are initiated by //-induced 
raising and lowering, respectively, of the {G(0\siS2)} relative to the fixed {G(l\siS2)}. The /i-induced change in the 
energy landscape cause the time-dependent population changes in responsible for the dissipative work of allosteric 
signaling. Population flux is represented as a ball rolling on the landscape. 




FIG. 1: Basin-limited portions of the free energy landscape {G(Si,/x)} of a two-component allosteric assembly (vertices), (a) 
System satisfying the 5 topological considerations (see text). Dominant pathways for activation (red arrows) and deactivation 
(blue arrows), (b) {G(Si,fi)} of the cardiac regulatory switch shown for desaturating Ca 2+ (fi — — 16.81fcsT), midpoint of 
activation (fi = — 13.31fcsT) and saturating Ca 2+ (/i = — 9.90/cbT). G(1\siS2) (Ca 2+ -bound species, black) are stationary; 
G(0\siS2, fi) (Ca 2+ -unbound species, red) depend linearly on /i = fcsTln[Ca 2+ ]. 

The topology of the basin-limited landscape {G(S^)} is constrained by five considerations, (i) The free eneries of all 
ligand-unbound species (0\siS2) depend linearly on fi: G(0\s\S2) = G°(0\sis 2 ) + n/x. n — 1 because cardiac troponin 
C has one regulatory binding site for Ca 2+ . The free energies of the ligand bound species G(l\siS2) are fixed (n = 0) 
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becuse the Ca 2+ binding site on troponin C is occupied, (it) The free energy changes along a closed path sum to zero, 
AGi^k + - ■ ■ + AGj < -i — 0, where AG^k = G(Si) — G(Sk)- Free energy conservation, along with microscopic detailed 
balance: TljiPi = TlijPj, produces macroscopically balanced transition probabilities, Tlik • • - Hji = Tlik ■ ■ ■ H_ji , where 
Pi = = P{(t — ► oo). (Hi) Long distance allosteric signalling is exponentially suppressed beyond the localization 
length of the low frequency modes that produce allostery [T7]. We assume that each low frequency mode is spatially 
limited to an assembly component and its nearest-neighbors. The assumption of nearest-ncighbor-limited influence, 
which parallels the Markov assumption, causes certain transition probabilities to be degenerate 

U(af\s . . . Sj-iSj+i . . . s m ) = U(of\8j-x8j +1 ) (2) 

(m = 2). The system-state Si — (sq\siS 2 ) has been rewritten as Si = (sj\sk^j) to emphasize the transition of 
component j — 0,1,2 against the fixed components k =/= j. Eq. ([2| implies, for example, that the reporter s 2 has 
no direct knowledge of whether ligand sq is bound to the sensor; rather, ligand binding is communicated to the 
reporter by allosteric change a\ in the sensor. Degenerate transition probabilities occur for activation/deactivation 
of the reporter crj, n(0\lo' 2 t ) = n(l\lcr 2 t ), and for ligand binding/release er*, II(<7*\si0) = H(a \sil). For these 
transitions, degeneracy is due to a degenerate free energy landscape and degenerate frictional coefficients upon which 
the Hj<—i depend. Particularly significant is degeneracy in the free energy change that governs switching in the 
reporter of, AG(0\lof ) = AG(l\lof ). Fig. [I] shows how nearest-neighbor-limited influence imposes parallelpipid 
geometry for of in the free energy landscape. Parallelpipid geometry is not imposed for of. (iv) Deactivation of the 
sensor while the reporter is active causes steric conflict, causing G(so\01) to be very high. The high energy states 
(so\01) are unvisited and can be ignored, (v) When ligand is bound, sensor activation is favorable AG(l\of 0) < 0; 
when ligand is not bound, sensor activation is unfavorable AG(0\of 0) > 0. 

Signaling error — We have used time-resolved and stopped-flow Forster resonance energy transfer (FRET) 
measurements [H] to parametrize the {G(Si)} of the cardiac regulatory switch: G(0\00, 0\10, 0\11) = 
{9.56, 12.52, 12.08} k B T+n; G(1\00, 1\10, 1\11) = {0, -2.89, -3.33} k B T; T = 15 C. All G{S l ) are relative to G(1\00), 
which is set to zero. Fig. [IJ> shows the parametrized landscape. Significantly, the free energy change that governs 
activation of the reporter is very small, AG(so\lc + ) = — 0.44fcsT. 

For the assembly equilibrated under the Ca 2+ chemical potential /i, the system-state probabilities Pi(fj) are Boltz- 
mann distributed 

Pi (fj,)=e-^ s ^/J2e- mSm) . (3) 
u 

Ca 2+ -induced signaling so — * 3 2 can be quantitated as the relative entropy change of the reporter S2 in response to a 
change in the Ca 2+ chemical potential \i [TH] 

V [{Ps 2 } \\{Ps 2 }] (M) = Y,Ps 2 ln fe/ftj , (4) 

{82} 

where p S2 = J2{ So Sl } Pso\sis 2 > s« = 0, 1 are marginalized system-state probablities. The reference distribution p* So — 
p S2 (/u — > — oo) is the equilibrium distribution in the absence of ligand. Systems with low signaling error mirror the 
input at the output (s2 = so) at all \x. There are two sources of signaling error: activity in the absence of bound 
ligand (constitutive activation) and inactivity when ligand is bound (incomplete activation). T> \p S2 \ addresses 
only error from incomplete activation. Both sources of signaling error £(/i) are treated by the conditional entropy 
H[s 2 \s ](fj) of s 2 given s [IS], 

£ = H[s 2 \s ](fi) = - P«o,s2 !og 2 [Ps ,s 2 /Ps ] , (5) 

{s ,s 2 } 

where p So , S2 = E{ sl }Ps \s 1S2 and Ps a = J2{ Sl ,s 2 } Ps \sis 2 - Signaling error, < £{p) < H[s 2 ] < 1, is a bounded 
logarithmic measure of the error in transmitting s — * s 2 for the cardiac regulatory assembly in equilibrium with an 
externally regulated fi. £q = H[s 2 \so]{— oo) = H[s 2 \sq = 0] is the error of transmitting sq = 0. £\ = ff[s2|so](oo) = 

H[s 2 \sq = 1] is the error of transmitting s = !• The mean signaling error £ = lim^^^ f* x dfJ,H[s 2 \so]/ J_ x d\x , 

which simplifies to £ = [£q + £\) /2, is a measure of system performance over all ligand concentrations. The fidelity 
of transmitting s$ = i is Ti = \ — £i. Overall system fidelity is T = J-qJ-i- 

The probability of being active p S2 =i an d signaling error £ were calculated as a function of fi from the parametrized 
{G(Si)} of the cardiac regulatory assembly. They are shown as solid lines in Fig. [2] The system has low constitutive 



FIG. 2: Signaling fidelity. Probability of being active p S2 = i Q (black), signaling error £ = H \s2\s0\ |5| (red), and minimum 
signaling error £ mm (green) vs. j3fj, — ln[Ca 2+ ]. Experimentally measured landscape (solid lines), artificial 1 ksT decrease in 
G(so\H) (dashed lines), and additional 1 ksT decrease in G(0\00) (dotted lines). 

activity, p S2= \{— 00) = 0.07, and incompletely activates (at high [Ca 2+ ]), p S2= i(oo) = 0.60. This generates reasonably 
low error at low [Ca 2+ ], £ = 0.37, but causes substantial error at high [Ca 2+ ], £ x = 0.97. The average signaling error 
is large, £ = 0.67. On paper, one can attempt to improve £ subject to the constraint that AG(0\lcr^~) = AG(l\l<rJ), 
abbreviated hence as AG. The dashed lines in Fig. [2]show p S2= i and £ after G(0\11) and G(l\ll) have been jointly 
decreased by 1 fceT. The modified landscape raises the maximum ligand-induced activity, p« 3 =i(oo) — 0-80 but 
also increases constituitive activity, p S2 =i(— 00) = 0.17. The desired reduction in £\ = 0.72 is offset by an increased 
£ = 0.66. In the re-designed system, £ = 0.69 actually increases. Recall that nearest-neighbor-limited influence 
does not constrain the free energy changes for the sensor, so G(0\00) and G(1\00) can be independently varied. The 
dotted lines in Fig. [2] show p S2 =i and £ after an additional adjustment to the landscape — lowering G(0\00) by 1 
ksT. The adjustment restores low constituitive activity, p S2= \{— 00) = 0.07 and retains the gains in maximal activity, 
p S2 =i(oo) = 0.80. £0 — 0.38, £\ — 0.72, and £ — 0.55. The modified landscape reduces average signaling error £ by 
18%. Below, we show that these gains in signaling fidelity are offset by a reduced rate of deactivation. 

What is the minimum possible signaling error £ mm for a given AG? An error-minimizing landscape with fixed AG 
populates (1\10) and (1\H) when Ca 2+ is bound and populates only (0\00) when Ca 2+ is unbound (i.e. £™ ln = 0). 
For all other species S{, = 0. The modified energy landscape of the cardiac regulatory assembly that produces 
£ min is /3G(0\00,0\10,0\11) = {-2.96, a, a - AG} + 15.40 + /3fi; and /3G(1\00, 1\10, 1\11) = {a + 2.89, 0, -AG}, 
where a — > 00. The error-minimizing landscape was subjected the same perturbations as above — 1 fc^T decrease in 
G(so\ll) and an additional 1 fcgT decrease in G(0\00). The green lines in Figj2]show £ mlD in the three landscapes. 
As required, £ > £ min for all p. Also, £^ in = and £f lin £ x . From ^ and ([5]), 

-/3AG 

= log 2 (1 + e-? AG ) + /?AG 1 + e ^ AG , (6) 

is the minimum attainable signaling error for a given AG > 0, and T niax = (1 - £ min )(l - £f m ) = 1 - £f m is the 
maximum attainable signaling fidelity. < :F max < 1 is bounded for AG > 0. 

Kinetics of resetting — The landscape of the assembly supports both activation and deactivation (resetting) stages 
of the signaling cycle. As shown in Fig. [lji, activation proceeds predominantly through the sequence — » o\ — > a^. 
Deactivation proceeds predominantly through the sequence — > — > cr^. The reporter deactivates before the 
sensor deactivates, causing the assembly to switch as a last-in-first-out (LIFO) stack. 

While a large AG increases maximal activation and lowers £\ , it also contributes to the free energy barrier (Fig [3] 
inset) that must be overcome (transiently) for deactivation to proceed. The reporter must gain —{T>b ~ T^a) = /3AG 
of information from the solvent before the sensor can deactivate. Here, Sa — (0\11), Sb — (0\10), and T>j = 
T>[S({Si} — Sj)\ \ {pi}] For a transition with free energy of activation AG^ = AG + 7^ and barrier crossing attempt 
frequency v, the absolute transition rate ks^A = vqUba satisfies the Arrhenius relation, k = vcxp [—/JAG"'"] = 
vexp [o^/fes] exp [— f3$~\ = v' exp [ — /3<5*] , where 5* is the enthalpy of activation, cfi is the entropy of activation, and 
v 1 is the entropy of activation-adjusted barrier crossing attempt frequency. The rate of information gain (AG > 0), 
-d T T> = ~(V b -V a )IIba is bounded: < -d T V < /3AGe" /3AG . See EPAPS Document No. [x] for proof. The upper 
limit is approached as v — > v$ and ^ — > 0. max(— d T T>) — e" 1 at /3AG = 1, and —d T T> — > as AG — > 00. Therefore, a 
large AG reduces signaling error (JsJ) but can critically slow deactivation. Indeed, temperature-dependent stopped-flow 
FRET measurements of the cardiac regulatory assembly [H] show that deactivation of the reporter (so\lo^~) is the 
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FIG. 3: Trade-off between the deactivation rate fc max and signaling fidelity _^" max j n ra te of computation I max by an allosteric 
assembly. 2T max = ^.max^-max j nge ^. £ ree ener gy surface for a network transition. 

rate limiting step in deactivation with i/' = 1.0 x 10 10 s" 1 and 5$ = 18.3 k B T (T = 15 C). Lowering AG by 2 k B T 

would slow (soVct^) b y = 117 -> L0 s_1 )- 

Physical limit The maximum rate of computation X max (in bits/sec) is the maximum rate of the information- 
gaining step in deactivation fc max = i/ Q exp(— (3 AG) times the maximum bit transmission fidelity jF max . We obtain 



with max(T max ) = 0.104^0 at AG = 1.12k B T. For a typical allosteric transition with wavenumber v = 10 cm -1 , 
Vq = cv = 3 X 10 11 s -1 , and ma,x(T max ) = 3.1 x 10 10 bits/s. This is about 1000 times less than the Hciscnbcrg limit on 
the rate of computation: 4k B T/h = 2.6 x 10 13 bits/s (T = 37 C) [13]. The Heisenbcrg limit ignores system resetting. 

Conclusions — Cells employ assemblies of allosteric proteins — molecular switches that bind diffusing ligands and 
communicate binding as altered activity — to regulate intracellular function. These protein complexes are nano-scale 
Brownian computers whose function depends on random collisions with solvent. Elucidating the physical properties 
that limit the rate of Brownian computation is essential for understanding the molecular systems-level design of 
biological signaling complexes and the pathophysiology of diseases that involve these complexes. The Ca 2+ -sensitive 
cardiac TnC-Tnl complex, an assembly of two allosteric proteins, is a paradigm for Brownian computation. The 
assembly functions as a digital logic buffer that transduces a binary input (the Ca 2+ bound/unbound status of the 
sensor TnC) to a binary output (the bi-metastable state of the reporter Tnl). 

We examined the cardiac TnC-Tnl complex using a non-equilibrium phenomenology of allostery in protein assem- 
blies. The phenomenology assumes limited spatial extent of protein-protein interactions within a protein assembly, 
just as residue-residue interactions have limited spatial extend within a protein |20[ I21j . Computation occurs as a 
Ca 2+ -induced perturbation of the the system's free energy landscape, causing a change in the population distribution 
among the set of coarse-grained system-states. Nearest-neighbor- limited influence, a spatial analog of the Markov 
assumption, constrains the topology of the system's free energy landscape and introduces degenerate transition prob- 
abilities. A single free energy change governs activation/deactivation of the reporter independent of whether Ca 2+ is 
bound to the sensor. Decreasing free energy drop enhances signaling fidelity but slows the rate of deactivation. The 
trade-off between fidelity and resetting speed physically limits the rate of information processing. This physical limit 
is an emergent constraint |22l 123] that arises from limited spatial extent of protein-protein interactions, a feature of 
the complex nature of the assembly. This limit is faced by assemblies that couple to a single chemical bath. Many 
assemblies appear to have addressed this limit by coupling to a second energy source. This energy source is usually 
inorganic phosphate derived from hydrolysis of nucleoside triphosphate. Examples include the family of G-protein 
coupled receptors, the NF-kB family of transcription factors, and the troponin-tropomyosin-regulated interaction of 
actin and myosin in striated muscle. 
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EPAPS 

Individual Markov transition probabilities ^Ij^i depend on the barrier crossing free energies AGj^ through the 
Arrhcnius relation pQ, 

f A^exp(-/3AG^ ; ), Vj^i 

with normalized barrier crossing attempt frequencies vJ^Zi = Vj^i/i>Q, vq = max({i/j<_i}). [3 — 1/k^T. The normal- 
ization constants Ni = Y]- Hj<—i ensure that ■ Ilj^j = 1, as required. All transition probabilities are bounded, 
< n„<_j < 1. The barrier crossing free energies are calcuated using 

7* + AG, V AG > 

7 t 



MJ U = \ ' 7 t ' AG < ( g ) 



where 7* > and AG = Gj — Gj. Furthermore, we set iVjzT^Ti = NjVjZ^i to satisfy microscopic balance, Ilj^j/IIj^j = 

exp(-/3G jV -i)- 

Substituting (|8| into Q we obtain 

_ _ f exp(-/3( 7 t + AG)), V AG > 
11,^ - NiUj^i I exp( _ / 3 7 t ); AG < ' (9) 

with the constraints > Ni,v^2i < 1 and 7^ > 0. Effective barrier crossing attempt frequencies (s _1 ) are defined 



PROOF OF THE EQUATION FOR THE BOUNDED RATE OF INFORMATION GAIN 

The relative entropy change for the transition po — > p from an initial distribution po to a second distribution p is 

£>[Po||p] = ^PoMpo/p]- 

u 

Systems will relax to the equilibrium density p* given by the Boltzmann distribution 

P*M = \e~^ s ^\ (10) 

where Z — exp —f3G(Si\ fi) is the partition function. Taking the system from state Sa = (0\11) with p = 5(S—Sa) 
to its equilibrium distribution p*, involves an information loss of T>a — 2?[p||p*] = — lnp^. Similarly taking the 
system from state Sb = (0\11) with p = 5(S — Sb) to its equilibrium distribution p*, involves an information loss 
of T>b — T>[p\\p*] — — Inps. The net gain of information for Sa —> Sb is — (£>b — £>a) = Mpb/pa]- Inserting (101, 
rearranging and canceling Z , the net gain of information is 

- (V B - V A ) = 0AG, (11) 

where AG = G(Sb) — G(S A ). When AG > 0, there is a gain of information. 

The rate of information gain for Sa — * Sb is defined, —d T V = —{T>b — T>a)TIb^a- Substituting (111 and ^ we 
obtain 

Finally, substituting for AG^ we obtain 

f)T>- RGN v / ex P(-^ + AG ))' V AG > 

-d T V - PGN a vb^a I exp( _ /37 t ); AG < ' 

But Na < 0, v~b<^a < 1 and 7* > 0. The rate of information gain is a bounded function (Fig. [4]), 

0< -d T V < /3AGexp(-/3AG), V AG > . , 

[3AG< -d T V, <0 AG<0 1 > 
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FIG. 4: The rate of information gain —d T T> for an allosteric transition with net free energy change AG. ( 12 1 (gray area 



The maximum possible rate of information gain occurs at the extremum d(—d T T>)/d/S.G = 0. Solving, we obtain 
max(— d T T>) — e _1 at (3AG = 1 as the maximum possible rate of information gain. There is no corresponding limit 
on the maximum rate of information loss. 

—d T T> should be interpreted as the rate of information gain per clock cycle. It is a dimensionless quantity because 
time has been normalized against the clock speed using r = ti/Q. 
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